Exact asymptotic expansions for the cylindrical Poisson-Boltzmann equation 
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The mathematical theory of integrable Painleve/Toda type systems sheds new light on the behav- 
ior of solutions to the Poisson-Boltzmann equation for the potential due to a long rod-like macroion. 
We investigate here the case of symmetric electrolytes together with that of 1:2 and 2:1 salts. Short 
and large scale features are analyzed, with a particular emphasis on the low salinity regime. Ana- 
lytical expansions are derived for several quantities relevant for polyelectrolytes theory, such as the 
Manning radius. In addition, accurate and practical expressions are worked out for the electrostatic 
potential, which improve upon previous work and cover the full range of radial distances. 
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I. INTRODUCTION 



Polyelectrolytes are polymer molecules bearing 
charged units. Within the associated relevant cylindrical 
geometry, the long range character of Coulombic inter- 
actions is responsible for the phenomenon of counterion 
condensation. In essence, the electric potential created 
by the charged polyion features a logarithmic depen- 
dence on radial distance. The competing entropy of 
confinement is of a similar functional form, which may 
result in a condensation of counterions onto the charged 
cylinder, when Coulombic interactions prevail, that is 
when the polyelectrolyte line charge exceeds a given 
threshold (see e.g. 0). Strongly charged Hnear poly- 
electrolytes thereby effectively reduce their line charge 
density. This was realized by Onsager in the 1960s 
and formalized subsequently " Q. It turns out that 
the mean-field Poisson-Boltzmann (PB) theory offers 
a valuable framework to discuss and analyze the phe- 
nomenon in detail 0, |1 i, II B IS ES IH iH III III 111 
Unfortunately, even within mean-field, exact results are 
scarce, and one has to resort to numerical resolution or 
propose approximations. 

This paper is dedicated to the derivation of exact re- 
sults within PB theory. We will be interested in the be- 
havior of a unique infinite charged cylinder of radius a 
and uniform line charge density e/6, where e > is the 
elementary charge and b may be viewed as an equivalent 
average spacing between charges along the polyion. This 
macroion is immersed in an infinite electrolyte contain- 
ing two species of microions (salt): their bulk densities 
ni and n2 far from the cylinder define the Debye length 
1/k through = AniBinizf + n2Z2). Here, Zit denotes 
the charge of species i and Ib = /3e^/e is the Bjerrum 
length, which involves the solvent dielectric permittiv- 



ity e and the inverse temperature (3. We will focus on 
the three situations where analytical progress is possible: 
1:1 electrolytes (or more generally symmetric z:z elec- 
trolytes), together with the 1:2 and 2:1 cases. We denote 
by 1:2 the situation where colons are monovalent (z+ = 1 
for a positively charged polyion) and counterions are di- 
valent (z_ — —2 or more generally, z_ = — 2z_|_). On 
the other hand, the notation 2:1 refers to divalent colons 
with monovalent counterions. To be specific, introducing 
the dimcnsionless radial distance f = Kr, the equations 
to be solved for the dimcnsionless electrostatic potential 
y — (}eip {ip being the electric potential) read 



(Py dy 

(Py dy 

df"^ f df 

d'^y 1 dy 

df"^ f df 



sinh(?/) 



1 {e-y - e-y) 



(1:1) (1.1a) 
(1:2) (1.1b) 
(2:1) (1.1c) 



The solution should satisfy the boundary conditions at 
a = K,a 



y ~ dy 

lim r — — 

f^a df 



-2C 



(1.2) 



and limf^oo y{f) — 0. We have introduced here the di- 
mcnsionless bare line charge ^ = that plays a key 
role below. It is understood that the cylinder is positively 
charged, without loss of generality |4^ . 

We shall not comment here the well documented limi- 
tations jjf the mean-field approximation (see e.g. [l^lT^ 
H?, 19, 20, "^Tl), but only summarize the main findings. In 
a solvent like water at room temperature, mean-field pro- 
vides an accurate description in the 1:1 case for all exist- 
ing polyelectrolytes. We also expect the 2:1 situation to 
be correctly described, while in the 1:2 case, micro-ionic 
correlations, discarded within mean-field, play an impor- 
tant role at large values of ^. Schematically, the coupling 
parameter F = z'^/^ \/C^B/(27ra) allows to appreciate the 
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importance of those correlations where z stands for 
counter-ions valency. When F < 2, mean-field holds. 
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The crux of our approach hes in the mapping between 
equations Hl.l|l and Painleve/Toda type equations (see 
below), that are in particular relevant to describe the 
spin-spin correlator in the two-dimensional Ising model, 
and have been the subject of deep mathematical investi- 
gations |23) 1^ • It seems that this body of work has 
not been fully transposed in the language of polyelec- 
trolyte physics, with the exception of the limit na 0^ 
(e.g. realized with a charged line -of vanishing thickness- 
at finite salt) where it allowed to show rigorously the exis- 
tence of the condensation 0, 0| and also provided exact 
results for the electric potential 0,^3. However, it will 
appear below that the limit na — > 0+ is approached log- 
arithmically slowly so that finite salt correction cannot 
be neglected in practice. It also turns that to the best of 
our knowledge, the case of 2:1 electrolytes has not been 
considered in the physics literature, whereas some results 
have been reported in the reverse 1:2 case, but below the 
condensation threshold and again in the limit a ^ |llj | . 
Our aim is to fill these gaps. Our formulas will prove 
extremely accurate when compared to the numerical so- 
lutions of Hl.lll . 

The paper is organized as follows. In section^ the for- 
mal solutions of Eqs. Hl.l|l will be given. As such, these 
relations are not useful in practice, and deriving "ready- 
to-use" and operational expressions will be the purpose of 
the subsequent analysis. Section ITTll will then be devoted 
to the derivation of asymptotic expansions for the elec- 
tric potential, and to the connection between short and 
large distance features. These results identify a change in 
the short distance behavior when increasing the charge 
density above a threshold ^c, what will be worked out. 
This is the fingerprint of counterion condensation, which 
turns out to be smoothed by salt (strictly speaking, the 
phenomenon is critical in the limit Ka —^ only) . Section 
llVl will analyze a few limiting cases of particular interest. 
Analytical expressions for several important quantities 
such as the Manning radius will be given and compared in 
section^to the numerical data obtained solving the non- 
linear PB equation. A similar analysis will be performed 
for the electric potential, which will assess the validity of 
the analytical results. Although particular attention will 
be paid to the low salinity regime (the results of sections 
IIIII tolV)relv on expansions that typically break down for 
Ka > 1), we will also report in section IVll some results 
concerning the effective charge of the polyion, valid for 
all values oi na. A preliminary account of parts of this 
work has been published in 25]. 

Before embarking on the analysis of the properties of 
Eqs. Hl.l|l . we note that for a macroion with a small 
reduced linear charge density ^cff, the problem may be 
linearized with solution 



yiin(r) = °,.. Ko{P). 
aKi[a) 



(1.3) 



Here Kq and Ki are the modified Bessel functions. Since 
the solution to the non-linear problem Hl.l|) vanishes 
when f — > oo we necessarily have that for large distances 



f, y{f) behaves as (|1.3ll . However in the non-linear case 
the constant feff is not the bare charge of the rod. This 
quantity is called the effective charge and depends on 
salinity conditions, bare charge, temperature etc. In the 
weak overlap approximation where the double-layers of 
different macroions are well separated (it is typically the 
case when their mutual distance is larger than 1/k), the 
(squared) effective charge governs the amplitude of the 
interaction free energy. 



II. FORMAL SOLUTION TO THE PROBLEM 

Poisson-Boltzmann equation Hl.lafl in the 1:1 case has 
been solved in Ref. in the context of the Painleve 
III theory, since e^/^ obeys a particular case of Painleve 
III equation (see e.g. |23| and appendix more details 
together with a relevant bibliography may be found in 
[231). More generally, Eqs. (|l.la|) . Ijl.lb|) and (|l.lc|) are 
the first equations of the hierarchy of cylindrical Toda 
equations. In Ref. a class of solutions to this Toda 
equations has been reported, hence the solutions to the 
Poisson-Boltzmann equations Hl.la|l . (|l.lb|l and (|l.lc|l . 
These solutions obey a boundary condition at f — > oo, 
given in terms of a constant A by 



for the 1:1 case, 



yii(f) ~ AXKoif) 



(2.1a) 



(2.1b) 



for both 1:2 and 2:1 cases. 

The solutions are expressed as determinants of certain 
operators. For 1:1 salts, we have 



yii(f) = 21ndet(l -f XKf) - 21ndet(l - \Kf) 
with Kf. an integral operator on with kernel 

Kf{u,v) = 



U + V 

In the 1:2 case the solution reads 

while in the 2:1 situation it is 
?;2i(r) -lndet(l + AKfV 
where 



Indet(l-Ai^f^) 
lndet(l + Ai^f^) 



iff = (c-c^)or + r-c)o^- 

(C'-i)of + (C-i)of 
(i-C)o.f + (i-C')of . 



(1) 



l(2) 



(1) 

(2) 



(2.2) 

(2.3) 

(2.4) 

(2.5) 

(2.6a) 
(2.6b) 
(2.6c) 



Here ^ 

on IR+ 



= e^"^^/^ and O-^-* and of''' are integral operators 
with kernel 



0f\u,v) 



-f[(l-C)«-(-(l-C')«-']/(2%/3) 

—C,u + V 



(2.7) 
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and 0^f\u,v) — 0\'''{u,v), the bar denoting complex 

conjugation. It can be shown [23| that det(l — AK^^^) = 

det(l- AiC^°^) (making use of a change of variable u — > 
thus from the solution to the case 1:2, yi2, one can 
obtain the solution for the case 2:1 as 2/21 = ~2/i2 with 
the change A ^ — A, as previously announced. 

To solve completely the problem we are interested in, 
we should impose the boundary condition (|1.2|l to ex- 
press A in terms of the bare charge ^. Notice that the 
constant A introduced above is closely related to the ef- 
fective charge ^cff- Indeed, 



U= 2dKi{&)X (1 
U= 3dKi{&)\ (1 
^eff= SdKi{d)X (2 



1) 
2) 
1) 



(2.8a) 
(2.8b) 
(2.8c) 



The determination of the effective charge in terms of the 
bare one is interesting per se, since it allows to use, at 
large distances, the linear theory expression, provided 
that the bare charge is replaced by the effective one. 
Writing down the boundary condition 111.2(1 gives the bare 
charge as a function of the effective one. In the limiting 
case a — > 0, this relation may be inverted since the associ- 
ated (so-called) connection problem was solved: knowing 
the large-distance behavior 1)2.1(1 of the solutions 1(2.2(1 . 
(12.4(1 and ((2.5(1 . the short-distance behavior for 7~ 
follows 0. The 1:1 and 1:2 cases for a = were studied 
in detail by Tracy and Widom in Ref. 'Tl'l. The 2:1 sit- 
uation, even for a = 0, was not explicitly considered in 
Ref. [ll| , although it can be obtained following the same 
lines as those exposed in [Tlj . 

The situation with a 7^ is analytically more diffi- 
cult, but we will show below that useful expressions may 
nevertheless be derived there, that rely on accurate ap- 
proximations for certain key quantities related to the pa- 
rameter A. These expressions allow to fully cover the 
regime of thin cylinders and/or low salt a < 1, whereas 
the discussion of results valid at arbitrary salt content is 
deferred to section IVTI 



Short distance behavior 



The general results of Ref. |2J| for the cylindrical 
BuUough-Dodd equation (which belongs to the Toda 
family and onto which the above asymmetric PB equa- 
tion can be mapped) allow to find the small distance 
asymptotics of ?/(f) given in Eq. ((2.4(1 for 1:2 electrolytes. 
Allowing for negative values of ^ will therefore also pro- 
vide the solution to the 2:1 case. Similarly, the 1:1 be- 
havior encoded in Eq. 1(2.2(1 follows from Theorem 3 of 
Ref. [22 and its corollaries (in particular those derived 
in section IV.I of [H). 

It turns that the parameter A plays an essential role 
here. Its position with respect to a threshold value Ac 
discriminates two different short scale behaviors. The 
threshold values read: 



1 

TT 



2\/3 7r 

73 

27r' 



(3.1a) 
(3.1b) 

(3.1c) 



We will show that the condition A < Ac is equivalent 
to ^ < £,c where ^c is some threshold charge. In the 
limit na — d —^ 0, coincides with the known Manning 
parameter beyond which counterion condensation sets in: 



^(1:1) 
*^Manning 



1:2) 



annmg 

^(2:1) 
'^Mannine 



= 1 
= 1/2 
= 1. 



(3.2) 
(3.3) 
(3.4) 



III. ASYMPTOTIC EXPANSIONS 



The results reported are obtained from the limit d 
but turn out to be reliable for a < 1 (see below). 



However, as soon as a 7^ 0, it will appear that signifi- 
cantly differs from ^Manning- For A < Ac, we have [a 113 



yii{f) ^ -2Alnf + 21nB-21n 
y^2{f) = -2Alnf + InS - In 1 
2/21 (f) -2Alnf + lnB - 21n 



1 - 



2 ~2-2A 



16(1-^)2 



1 - 



12(2^1-1)2 

~2-2A^ 



24(1-^)2 



+ 0(f2+2^) 

+ 0(7=2+2^) 

+ 0(f2+4^) 



(1:1) 
(1:2) 
(2:1) 



(3.5a) 
(3.5b) 
(3.5c) 
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where A is a function of A, see Eqs. (|3.11ll and H3.13|l below, and 



B 



B = 



B = 



l-2A^ 



q3ylr,2A ^13/^1 3 I 



(1:1) 
(1:2) 
(2:1) 



(3.6a) 
(3.6b) 
(3.6c) 



Here F is the Euler gamma function. 

On the other hand, for A > Ac (or equivalently ^ > ^c), 
the short distance behavior is characterized by a param- 
eter /i > and we have 



4/. 



-2^1og(f) - 2^C(^^^^1 (3.7a) 



3\/3/^ 
V2f 
3V3 



-3/ilog(f)-3AiC(i^2)j (3 7b) 
(3.7c) 



|^log(f)-^MC(^^^) 



For the sake of notational simphcity, we do not exphcitly 
mention that ji depends on the situation 1:1, 1:2 or 2:1 
considered. The constants C appearing above are defined 
as 



^(1:2) 
^(2:1) 



7- 31og2 ~ -1.502 

1 3 
7 — - log 2 — - log 3 c 



7 - log 2 - - log 3 



(3.8) 

-1.301 (3.9) 
-1.763 (3.10) 



where 7 ~ 0.5772... is the Euler constant. We empha- 
size that expressions (|3.7|l do not have the same status 
of mathematical rigor as their counterparts 1)3. 5|l valid 
for A < Ac, as appears in Appendix IbI The constants 
C in (|3.7|) arise from the linearization of a /i-dependent 
function, which is justified from a practical point of view 
since for physically relevant situations, fi is small enough 
to allow for the corresponding Taylor expansions. For 
the exact expressions see Appendix 151 

We note here that the relation between surface poten- 
tial and surface charge follows from the expressions given 
in the present section. At this point, the electric potential 
depends on a parameter A for A < Ac or on a parameter 
/i for A > Ac. These two parameters are related to the 
bare reduced charge f of the cylindrical polyion through 
the boundary condition H1.2|l . Before clarifying this rela- 
tion, we precise the connection between small f and large 
f behavior. 



B. Connecting short and large scale features 

For 1:1 electrolytes, the connection problem was solved 
in H 112, with the result 



1 . /ttAA 
\ — — sin for 

A = — cosh(7r/i) for 



A < vr^i (or i < 6) (3.11) 
A > TT-i (or e > ^c)(3.12) 



Since A is related to the effective charge with governs 
the far field behavior [see H2.8|) ]. the above expressions 
realize an interesting connection between small and large 
f features. 

For asymmetric electrolytes and A < Ac 



A= -- 



1 



3 . fl 

2^™ 2 



13 . / 1 

A— - + — arcsm — - 
4 27r V 2 



A 
2A^ 
3A 
2A7 



(1 : 2) (3.13a) 
(2 : 1) (3.13b) 



In these two expressions, the values of Ac differ [A^^'"'^'' = 
3X'^'^^ see ()3.1|l ] in such a way that expressing A as a 

(1-2) 

function of Ac in both cases would provide the same 
expression, but a change in sign in A and in A. This 
reflects the original symmetry of 1:2 and 2:1 situations, 
when ^ is allowed to change sign. This symmetry is bro- 
ken here by the choice ^ > in both cases, which illus- 
trates the practical difference between them. 

From a more formal point of view, the expres- 
sions H3.11|l and H3.13|l hold even for complex parameters, 
provided by 



A ^ (-(X3,-A(i^i)]U[A(i^i),oo) (1 

A ^ (-(X3,-A(2^i)]U[A(i^2),oo) (1 

A ^ (-(^,-A(1^2)]u[^(2:l)^^) (2 

Alternatively, for A > Ac, we have 



-^ + 1 = 2 cosh(7r/i) 
Ac 



3A 



2 cosh(7r/i) 



(1:2) 
(2:1) 



1) (3.14a) 

2) (3.14b) 
1) (3.14c) 



(3.15a) 
(3.15b) 
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C. Applying the boundary condition at polyion 
contact 

In order to have a closed problem, we need to impose 
the boundary condition H1.2|l at f = a which relates A to 
^ (and therefore A to ^ from the connection formulae) if 
A < Ac and similarly relates /i to ^ if A > Ac. 

A straightforward computation gives in the 1:1 situa- 
tion 



{2-2A){Ka) 



2-2A 



16(1 - A)^B^ - (ko 
1) tan [2fi log(Ka/8) + 2^7] 



2^ (3-16) 
2/i (3.17) 



where it is understood that the first line holds for A < Ac 
while the second is relevant for A > Ac, in which case fi 
is the smallest positive root of the equation. Here again, 
7 denotes the Euler constant. The dependence of ^ and 
A upon ^ is shown in Figures ^ ^nd [3 




FIG. 1: Short scale exponent A (relevant for A < Ac) as 
a function of bare charge for a 1:1 salt. Four salinities are 
displayed: from left to right, a = 10"^ 10"^, lO"'^ and 0+, 
the latter coinciding with the first bissectrix. The values of 
^ where A = 1 correspond to an end point which defines the 
threshold bare charge ^c. 



In addition for the 1:2 situation, we get 
(l-2A)(Ka)2-4^B2 



A 



12(2/1- 1)2 -B2(Ka)2-4^ 
(2^-1) tan 3fi\og{Ka) + 3nC^^'-^'> 

while for the 2:1 case 

2(1- A)(Ka)2-2A5 



24(1 - yl)2 - (ko) 



2-2A 



B 



(^ - 1) tan 



-Mlog(«a) + -/.C(2^i) 



(3.18) 
= 3/i (3.19) 

(3.20) 
= |m(3.21) 



The values of B are given in H3.6|) . 

At this point, we emphasize that the above expressions 
are not mathematically exact, since use was made of the 
asymptotic expansions given in section Fill Al to compute 
dy I df in 1)1. 2|l. They nevertheless become asympt otically 




FIG. 2: Short scale exponent jj, relevant for A > Ac = I/tt 
versus reduced bare charge (1:1 salt). The three curves cor- 
respond from left to right to a = 10"^, 10"^ and 10"^. The 
limit 5^0 gives fi — 0. The salt dependent threshold val- 
ues ^c are shown by the arrows, which correspond to the end 
points where A = 1 in Figure Q The dashed curves corre- 
spond to approximation I4.13at . 



exact for na = d ^ [45j| . and will be shown to provide 
excellent results for the full range of thin cylinders a < 1 . 
Finally, it is important to remember that 



A *■ '^Manning 



for A 
for A 



K 
Aj 



(or e 
(or^ 



C) (3.22) 
O- (3.23) 



In the vicinity of the threshold ^c 
^-^Manning and ^l scalc like 
from Figures and [5] 



it can be checked that 



^c, as may be expected 



D. What are the associated threshold potentials 
and charges ? 

To understand the change of behavior of the electro- 
static potential below or above Ac, it is instructive to 
compute the associated threshold charge £,c- There are 
several ways to perform such a computation, either from 
below taking the limit A ^ A~ and considering the re- 
lations between A and ^, or from above (A A+) ma- 
nipulating ij,. By construction (see Appendix ^ for the 
details of the analytic continuation method used), these 
two routes provide the same short distance potentials, 
that read for the three different electrolytes 
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yii(^) 



yi2ir) = 



2/21 (r) 



-21og( - ) -21og 



-log(f) -C^i^i) 



-log(^) - log -log(f)-C(i^2) 



-21og(-^) -21og -log(r~)-C(2^i) 



(3.24a) 
(3.24b) 
(3.24c) 



r 



Limiting ourselves to the region a < 1 (as for all the 
results reported in this section ^Oj) ensures here that 
the previous expressions are real. Equations ( |3.24a|l and 
H3.24b|) were already given in |^ 0, [S^, l24| and are 
repeated here for the sake of completeness (expression 
(I3.24b [l appeared under number (3.2) with misprints in 
Ref. jllf. where the notation 21 used there corresponds 
to our 1:2 case [i^'l. 

From Eqs. H3.24|l . it is straightforward to compute the 
integrated charge q{r) in a cylinder of radius r which 
reads q{f) — —{f/2)dy/dr from Gauss' law. Evaluating 
this expression at f = a, we recover ()1.2|l . hence the value 
of ^ associated to Ac 



^(1:2) 



= 1 



1 



log(Ka) +C(i^i) 
1 1 



= 1 



2 21og(Ka) +2C(i^2) 
1 



log(Ka) +C(2^i) 



(3.25a) 
(3.25b) 
(3.25c) 



where the constants C are given in H3.8|l . We will see 
in section IIV Al that these threshold values, which dis- 
criminate A < Ac behavior from that at A > Ac, can be 
associated with the phenomenon of counterion condensa- 
tion. 

Expressions H3.25|l correspond to the leading order be- 
havior of f c and may be improved from the accurate ex- 
pressions obtained in |22. l24l|. In the 1:1 case, inclusion 
of the so far omitted dominant correction leads to 



29 



0(f5 log2(r~)), (3.26) 



with 



= -iog(f) -c^i^i). 



(3.27) 



Equation H3.24a|) corresponds to truncating the right 
hand side of (|3.26|l after the first term. The integrated 
charge q{r) associated to (|3.26|l reads, including only the 
dominant correction to the q{r) leading to H3.25a|l 



q{f) 



1 



1 



16 



+ 0{f^n^if)). (3.28) 



From q{d) = ^, we get that the error made in (|3.25a|l — 
which corresponds to the first two terms on the right 



hand side of (|3.28|) evaluated at f = S, — is of order 
^■^(loga)^. More quantitatively, the terms neglected in 
(|3.25a|l are below 2 x 10^'^ and are therefore irrelevant 
given that the term in 1/ log d in H3.25a|l may be of order 
0.1 or larger under reasonable salt conditions. A similar 
analysis could be performed to improve over expressions 
(l3:2tbll and JiO^ . 

Finally, note that exactly at ^ = ^c, effective charges 
take a particularly simple form [making use of Eqs. H2.8|l 
and H3.1|) ] since then A = Ac. 



IV. A FEW LIMITING CASES 

A. a = and arbitrary ^ 

It is instructive to discuss first the case a = (at fixed 
K > 0) where the analytical solutions take simple forms 
and A can be found explicitly. From Eqs. (|3.16|) . (|3.18|) 
or H3.20|l we find A = ^ and 



A = 

A = 
A = 



2^/37^ 
1 

2V37r 



2 sin 



2 sin 



1 . 






— sm 






TT 


2 




/2tt 


(- 










4 


/2tt 


(- 










4 



- 1 
+ 1 



(1 : 1) (4.1a) 
(1 : 2) (4.1b) 

(2 : 1) (4.1c) 



From these results we obtain the effective charge of the 
rod. For a = 0, we simply have [see Eq. H2.8|l ] ^off = 2A 
(1:1) and Coff = 3A (1:2) or (2:1), thus 



^eff 
Ccff 
^eff 



2 . TT^ 

— sm — 

TT 2 



V3 
2tt 
V3 
2n 





'2n 








2 sin 




i)] 


- 1 




Y 








'271 








2 sin 








+ 1 




Y 









(1 : l)(4.2a) 
(1 : 2)(4.2b) 

(2:1) (4.2c) 



Eqs. I4.2a|l and H4.2b|) may be found in 9] and j^. 
These formulas are valid only if the conditions H3.14|l on 
A are satisfied. Given that, for a = 0, the threshold 
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charges (|3.25|l read 

^(1:1) = 
^(1:2) ^ 



(1:1) 
anning 



^(1:2) 
^Manning 



(2:1) 



anning 

expressions H4.2|l hold provided 

e<i (1 

^<l/2 (1 
C<1 (2 



= 1 

_ 1 

^ 2 

= 1 



1) 
2) 
1) 



(4.3a) 
(4.3b) 
(4.3c) 



(4.4a) 
(4.4b) 
(4.4c) 



or, in the most general case ahowing for negative values 
of the bare charge: 



- 1< e < 1 
-1 < e < 1/2 
-1/2 < e<l 



(1 
(1 
(2 



1) 
2) 
1) 



(4.5a) 
(4.5b) 
(4.5c) 



We recover here the Onsager-Manning-Oosawa criterion 
for counterion condensation 0, |3j S H^- 
Strictly speaking, if ^ does not satisfy H4.4|l there is no 
physical solution of the Poisson-Boltzmann equation for 
a = 0. As explained in if the bare charge ^ is above 
the Manning-Oosawa threshold given in Eqs. H4.4|) . the 
Boltzmann factor of the interaction between the rod and 
a counterion is not integrable at short distances, and thus 
would lead to a collapse between the macroion and the 
counterions. The situation is similar to the one in the the- 
ory of two-dimensional Coulomb systems [2^ , with 
logarithmic interaction, where a system of point particles 
is stable only if the charge of the particles is small enough 
so that the Boltzmann factor of the interaction between 
two oppositely charged particles is integrable. Other- 
wise only a system with hard-core particles (or with other 
short-distance regularization of the Coulomb potential) 
is stable. 

For the 1:2 electrolyte, the value ^ = CMarrning = 1/2 
corresponds to Manning-Oosawa threshold |l,y|- Beyond 
the threshold (^ > 1/2), the effective charge ^cff attains 
its saturation value, — \/3/{2tt) ~ 0.275. One may 
check that in the linear regime (weak charges, ^ <C 1), ■^off 
and ^ coincide, as they should. For the 2:1 electrolyte, 

1 



(2:1) 



the Manning-Oosawa threshold is for ^ = ^Manning 
and beyond this threshold the effective saturated charge 
is given by — 3\/3/(27r); it is three times larger than 
in the 1:2 case. In the symmetric 1:1 electrolyte case. 



j(l:l) 



= 1 



the Manning-Oosawa threshold is for ^ — ^Manning 
and the corresponding effective saturated charge is ^^|* = 
2/7T ~ 0.636. We note that this result is in harmony with 
exact bounds derived by Odijk in [s^ : 0.59^^^* < 0.67. 

This analysis clarifies the relationship between two dif- 
ferent notions, the Manning-Oosawa threshold for coun- 
terion condensation and the notion of effective charge. 
The key point here is that the Manning-Oosawa thresh- 
old ^Manning for the bare charge is different from the 



value of the saturated effective charge ^^|*. The fact 
that ^^1* < ^Manning IS of course a non-linear effect 
which means that even accounting correctly for counte- 
rion condensation, the remaining layer of "uncondensed" 
or "free" microions cannot be treated within a linearized 
(Debye-Hiickel-like) theory. This exact result contrasts 
with common belief in the field, which possibly takes its 
roots in the work of Manning 2J . 
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FIG. 3; Effective vs. bare line charge for an infinite charged 
rod (1:2 case). The analytical result (I4.2bl valid in the low 
salt or thin rod limit kq 0, shown by the thick continu- 
ous curve is compared to the numerical solution of Poisson- 
Boltzmann theory for na = 10~® and Ka = 10""'. 
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FIG. 4: Same as Fig. Elfor the 2:1 electrolyte. Inset: Ccff/C 
for small values of ^. Notice the initial overshooting effect 

Ccff > i- 



In Figs. 13 and 01 we compare the result of Eqs. I)4.2b(l 
and H4.2cl) for a 1:2 and 2:1 electrolyte, respectively, with 
numerical data, for two low values of na. All numeri- 
cal PB data have been obtained following the method 
discussed in [sj- One may conclude that the limiting 
behavior ko ^ is reached very slowly, in fact, logarith- 
mically, as we will show below. 
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It is interesting to notice in Fig. 0] that in the 2:1 
electrolyte we have the overshooting effect, previously 
reported in Ref. , where the effective charge becomes 
larger than the bare charge for intermediate values of the 
latter. Indeed for |^| ^ 1, from Eqs. H4.2|l . we have for a 
2:1 electrolyte 

eeff = C+^-^+0(^^) (2:1) (4.6) 

and for a 1:2 electrolyte 

<2 2n^e 



3%/3 27 



+ 0(e) (1:2) (4.7) 



The first deviation of the effective charge from linear be- 
havior is positive in the 2:1 case (overshooting effect) and 
negative in the 1:2 case (no overshooting). For a 1:1 elec- 
trolyte, 



(4.8) 



The first deviation of the effective charge from the bare 
one is of order (as in the case /ta » 1 |22|) and negative 
(thus no overshooting). There is no term in as opposed 
to the charge asymmetric electrolytes 1:2 and 2:1. The 
physical origin of the overshooting effect lies in the fact 
that in the 2:1 case, the divalent colons are expelled fur- 
ther away in the double layer than the monovalent ones in 
the 1:1 situation, which results in a stronger electrostatic 
potential at large distances. Similarly, the screening for 
1:2 electrolytes is more efhcient since divalent countcri- 
ons are available, which not only results in the fact that 
Ceff < but also in the lower value of . We therefore 
speculate that the overshooting effect is generic in a z:z' 



electrolyte, with colons of valency z 
counterions (z'). 



larger than that of 



example for the 1:1 case, we obtain, if A < 1, 



,-Vll/2 ^ 



1 - 



16V(1-^)^ 52(1 -hA)2 



+0(f2(2+2A)^ 



(4.9) 



Notice the expected symmetry in the two terms when 
one changes A — > —A (in which, from the definition of 
B, gives B — > B^^). Taking the limit yl — > we obtain 



yii(^) 



-2 A 



7 + In - - 1 



(4.10) 

We recognize the small-r expansion of the Bessel function 
Koir). On the other hand l(TTT|l says that for A —> 0, 
A ~ 2X ^ ^cff/(aA'i(a)). We conclude that the small-r 
expansion l|4.10|l is the small-r expansion of the linear 
solution 11. 3|) and that, clearly, imposing the boundary 
condition (|1.2|l at r = a will yield ^ = ^eff- 

A similar conclusion is reached for the asymmetric 
cases. For the 2:1 case the small-r expansion, for A < 1, 
reads 



3-2/21/2 ^ 



1 - 



-2+4^ 



24(1 -yl)2 2452(2^1+1)2 



+ 0(f 2(2+4^)) 



(4.11) 



while the one for j/12 can be obtained from the preceding 
expression through the substitution A ^ —A and B 
B~^. It is direct to verify that in the limit A ^ Q, the 
asymptotics H4.11|l yield again H4.10|l . that is the small-r 
asymptotics of ?/i2 and y2i are those of the Bessel function 
Koif), solution of the linear problem. 



B. C ^ and arbitrary a < 1 

In the limit where ^ 0, the solution H1.3|l becomes 
exact for all distances, and one has £,cs/£, 1, as was 
simply checked for a = in section HV Al At finite salt 
however, the requirement ^cff/C ~^ 1 when ^ — *■ provides 
quite a non trivial benchmark for our analytical expres- 
sions. It turns that the results of sections llll Bl and 1111 CI 
are not sufficient -at finite a— to show fcff/C ~^ 1- This is 
because in Eqs. H3.5|l we neglected terms of order ^2+2-^ 
in the 1:1 and 1:2 cases and terms of order r2+4A j^-^ ^j^g 
2:1 case. As ^ ^ 0, we have A — 0, thus these terms be- 
come of order r2 and are of the same order of the last term 
in the logarithms of Eqs. (|3.5|l . Then, to properly com- 
pute the short-distance asymptotics as f ^ 0, we need 
the next order terms in Eqs. (|3.5|) . These can be easily 
obtained by replacing the small-r asymptotics H3.5|l into 
the Poisson-Boltzmann differential equation (|l.l|l . For 



C. C = ^Manning and arbitrary a < 1 

It is in general not possible to find analytically the so- 
lution to the transcendental equations satisfied by the 
parameter /i, which quantifies the behavior at short dis- 
tances above the salt dependent threshold, ie for ^ > ^c- 
However, if ^ = ^Manning (which is always beyond ^c), the 
relevant roots of Eqs. (|?T7|) . l|TT^ and read 



,(1,1) - _ 



,(2,1) 



TT 


4 [log a + 




TT 




6[loga 4 


C(i^2)] 


TT 




3 [log a + 


C(2:l)]- 



(4.12a) 
(4.12b) 
(4.12c) 



These expressions will be useful in section Ivl 
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D. 



^ > ^Manning and arbitrary a < 1 



The limit of large enough ^ allows to derive a use- 
ful approximation for /u, that will turn important from 
a practical point of view. Coming back to Eqs. H3.17|l . 
(|3.19|) and H3.21|l . we note that for diverging ^, the argu- 
ments of the tangent functions should be close to — tt, so 
that the tangent functions vanish, to be compatible with 
a finite value of /i. Expanding then the tangent to first 
order, we obtain 



(1,1) 



,(1,2) 



(2,1) 



TT 


2 [log 




i/(e-i)] 




TT 




3 [log 


;a-fC(i^2) _ 


1/(2^-1)] 








3 [log 







(4.13a) 
(4.13b) 
(4.13c) 



The domain of validity of these relations is estimated to 
be 1^ > ^Manning + C(l/loga), as may be observed in Fig. 
|21in the 1:1 case. In the vicinity of ^Manning where (I4.13f) 
fails, one should resort to expressions H4.12|l . We finally 
note that when ^ diverges, /i saturates to a finite value 
[twice those reported in (|4.12|l ] 



(1.1) _ 

Msat — 

(1.2) _ 
Msat — 



(2,1) _ 
Msat — 



2[log&-HC(i^i)] 

TT 

"3[log&-hC(i^2)] 
27r 

"3[loga-f-C(2:i)]- 



(4.14a) 
(4.14b) 
(4.14c) 



From Eq. H3.12|l and (|3.15|l we conclude that A does also 
saturate: a corollary is that effective charges -directly 
related to A throug h (j2.8l) - saturate, a generic feature of 
mean-field theories \32 



RESULTS AT FINITE SALT AND 
DISCUSSION 



The results reported in sections IIIII and IIVI are based 
on asymptotic expansions in the limit a — > 0. This limit 
is approached logarithmically slowly, so that finite d cor- 
rections should always be important in practice. We now 
therefore address the question of the reliability of our 
expressions at finite salt, by confronting them with the 
numerical solution of the non-linear PB equation Hl.l|l . 



A. Counterion condensation and Manning radius 

One of the most interesting features emerging here is 
that it is possible to generalize the notion of counterion 
condensation to finite salt systems. It is important how- 
ever to emphasize here that condensation is not an all 
or nothing process, that would occur precisely at ^c- In- 
stead, it is a gradually built phenomenon. This is partic- 
ularly true at finite salt density but already in the limit 



of vanishing salt Ka ^ 0, it is noteworthy that non linear 
effects play an important role for C < ^c, as revealed by 
the fact that ^eff/C 1- 

We concentrate here on a 1:1 salt, but the same anal- 
ysis also holds for 1:2 and 2:1 salts. The value dis- 
criminates between the short distance behaviors (|3.5a|) 
for £, < £,c and (|3.7a|) for ^ > ^c- As may be ob- 
served in Fig. n A is close to ^ except in the vicinity 
of the threshold charge ^c, so that to dominant order, 
yii behaves as — 2^1ogf, which corresponds to the po- 
tential of a bare cylinder with reduced charge ^, un- 
affected by screening. At = ^c, the potential de- 
velops an additional log log r term (see section IIII D|) 
while for £, > £,c, the behavior stemming from (|3.7a|) is 
more complex. For high enough ^ (more precisely, for 

C > ■^Manning + \O{l/l0ga)\ > ^Manning (which is itSClf 

larger than ^c), one may resort to approximation (|4.13() 
and after some manipulations, rewrite (|3.7a|) as 



-!/ll/2 



~ 2 



log 



e-1 



(5.1) 



Here, we have used the fact that 



2^ (logf + C(i^i)) = -7r + 2/ilog(0 



2/x 



so that Eq. holds for log(r/a) < 0(l/u].. Eq. 



(5.2) 

EU 

had already been derived by Ramanathan [T] (see also 
appendix A of 'S?!) and implies that to dominant or- 
der, the potential behaves like —2 logr at short distances, 
which corresponds to the bare potential of a polyion with 
C — ^Manning = 1, and is the fingerprint of counterion con- 
densation. We also note that for high ^, Eq. H5.1|l yields 
a total concentration of ions close to the rod that is pro- 
portional to the square of the surface charge density, as 
in the planar case [33 ■ To complement the above results 
that assume f > CManmng+ |C'(1/ loga)|, we also mention 
that for ^ = ^Manning (iG abovc the threshold), the results 
of section ITV CI yield 



-m/2 ^s^- _Ifloga + C(i^i; 

TT V 



(5.3) 



The dominant behavior for yn is therefore again —2 log r. 
Finally, decreasing further ^ to investigate the regime 
where it is close to the threshold (but still larger than 
^c), one may take advantage that /i vanishes at ^ = to 
get Eqs. 



-yil/2 



-2/21/2 



(logr + C(i^i)) (5.4a) 



= -_ (logf + C(i^2)) (5.4b) 



V6 



logf + C(2^^' 



(5.4c) 



The dominant behavior at short scales is again y ~ 
— 2^Manninglogr, which is the unscreened potential cre- 
ated by a rod of reduced charge ^Manning- This is quite 
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remarkable given < '^Manning since it holds in partic- 
ular for a cylinder with ^ verifying (^^ <)£. < ^Manning- 
Such a remark is nevertheless quite misleading since at 
finite values of Ka, the requirement r > a does not allow 
to take the limit r ^ where the different leading and 
sub-leading contributions to the potential can be identi- 
fied: "corrections" to the "leading" term are important. 
In other words, the full expression is required to approxi- 
mate y in H5.4|l and the "dominant" term — 2^Manning logf 
is in practice a bad approximation. 

Our analytical expressions identify a mathematical 
change in the behavior of the electric potential at ^ = ^c, 
that may be considered as being associated with coun- 
terion condensation. Such a terminology may however 
be confusing since at finite a, no singularity signals the 
crossing of the threshold value This restriction should 
be borne in mind. Of particular interest is the condensate 
structure, completely encoded in Eqs. 113. 7|l . and in par- 
ticular the condensate thickness [ij, uM ■ definition 
is necessarily arbitrary. The so-called Manning radius is 
often considered when it comes to quantifying the con- 
densate size. It is defined as the distance where the 
integrated charge g(r) = —{r/2)dy/dr equals ^Manning- 
Since we have seen that S,c < ^Manning? it is clear that such 
a criterion in inoperant for < C < ^Manning, which con- 
stitutes quite a deficiency, but the corresponding value 
nevertheless exhibits remarkable features within its 
domain of definition ^ > ^Manning- We first compute 
q{r), that follows directly from Eqs. H3.7|) 



Km -. 



g(2:l)(^) 



1 + 2/i tan 
1 3^* 



tan 



3n 

1 -I- tan 



-2/ilog 



-3^ log 



„(i:i) 



„(1:2) 



3/i 



(5.5a) 
(5.5b) 
.(5.5c) 



with 



(1:1) 



,(1:2) 



,(2:1) 



exp 



exp 



exp 



C(l:l) 


-- 


(5.6a) 


^(1:2) 




(5.6b) 


^(2:1) 


--) 

3^7 


(5.6c) 



In these expressions, one observes that q{r) = ^Manning 
for r = Tm, which means that the r„i appearing in H5.5|l 
are by definition Manning radii. When ^ = ^Manning, 
then by definition the Manning radius should coincide 
with the polyion radius a. This is indeed the case here, 
as a consequence of Eqs. (|4.12|l . For large enough ^, we 
obtain an expression in closed form from approximation 



,(1:1) 



,(1:2) 




2 2{i-l) 
In the cell model, the previous definition 

Q'(^m) — ^Manning 



(5.7a) 
(5.7b) 
(5.7c) 

(5.8) 



offers a geometric construction to compute rj„: q{r) plot- 
ted as a function of log r displays an inflection point at 
r = r,„ (see Ref. [H where it was also shown that this 
criterion could be extended beyond mean-field and used 
in Molecular Dynamics simulations to define a fraction of 
condensed ions). From the functional form of Eqs. H5.5|) 
-with a shifted tangent as a function of logr, i.e. the 
same form as in the salt-free cell model case- we see that 
a similar inflection point criterion holds here. This will 
be illustrated further in section IV CI 

For a 1:1 electrolyte, (|5.7a() is precisely the result ob- 
tained in Ref. [l5l |. The present work therefore pre- 
cises the domain of validity of this expression, and offers 
with H5.6a|) supplemented with (|3.17|) a better approx- 
imation (see Figure [5|. By construction, the accuracy 
of our expressions improve upon decreasing a. To test 
the worst cases, we therefore chose a relatively "high" 
value d — 0.01 in Fig. [S] where the agreement is seen 
to be very good, while (|5.7a() fails when ^ is too close 
to ^Manning = 1- Fig. iUustratcs the effect of increas- 
ing salinity and displays the results associated to three 
values of a. It is observed that at high ^, the analytical 
predictions deviate all the more from the numerical re- 
sults computed for PB theory as d is increased. For the 
highest value a = 0.3 of Fig. El the error made when 
^ ^ oo is on the order of 30%. It also appears that 
even at relatively high a, the analytical prediction is re- 
liable for low ^ [bearing in mind that ^ must be larger 
than ^Manning to allow for the deflnition of rm through 
Eq. (|5.8|) ]. For completeness, we present results for the 
1:2 case as well in Figure [7| changing a at flxed charge 
^ = 0.55 which is just above the value ^Manning = 1/2. 
The results at saturation are also reported (dashed line 
and squares). One needs to push d beyond 0.1 to no- 
tice a difference between (|5.6b|) and the numerical PB 
result (see the inset). Figure [T] illustrates how finite salt 
effects influence the scaling exponent a in the relation 
(wa)". If C > ^Manning + jC' (1/ log a) | , that is at 
high enough ^ or low enough a, we have a = 1/2 as dic- 
tated by Eqs. H5.7|) and evidenced by the dotted line in 
the main graph. If on the other hand one sits in the re- 
gion ^Manning < <^ < CManning+|C(l/loga)|, the exponent 
a changes continuously, and increasing d at fixed ^ drives 
the system toward the regime a = 1, a feature -visible 
for the data at ^ = 0.55 in Fig. [T]- which simply refiects 
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the fact that — > a. Finally, we emphasize that the 
behavior in all three 1:1, 1:2 and 2:1 situations are quali- 
tatively very similar so that the conclusions reached and 
phenomena observed are transferable from one situation 
to another. 




□ 2 4 pe 8 10 

I ' ^3 

1.0 2.0 



FIG. 5: Manning radius as defined from Eq. 115. 8|l versus 
reduced bare charge 4, for a 1:1 salt with na — 10~^ (a is 
the polyion radius). The continuous curve shows the results 
obtained from the numerical solution of PB equation, while 
the circles stand for expression I5.6all where ^ is the smallest 
positive root of 13.171 (see Fig. |21for a plot of as a function 
of bare charge). The squares indicate the prediction 15.7all 
(see the dashed line in Fig. |21 for the associated fS) which 
becomes asymptotically correct for large ^. The inset displays 
the corresponding charge regime, where rm/a, saturates at 
high 5 to the value 2^2/aexp(— 7/2) indicated by the dashed 
horizontal line. 



1.0 



0.5 



Kfl = 0.3 




A 


A 


A, 


A 










Ka = 10~' 








Ka = 10~' 








—f\ 


A 





1 2 3 4 5 

FIG. 6: The Manning radius r™ versus the reduced bare 
charge ^ for a 2:1 salt, for different values of na. The con- 
tinuous curve shows the prediction of Eq. 15.6cll while the 
triangles correspond to the numerical solution of PB theory. 

The inflection point feature stemming from (|5.8|) has 
the merit to unify the present infinite dilution/finite salt 
phenomenology with the finite density /vanishing salt sit- 




Ka 



FIG. 7: Manning radius as a function of salinity a for two 
different charges and a 1:2 salt. The symbols show the pre- 
diction of Eq. 115. 6bt (triangles for ^ = 0.55 and squares in the 
saturation hmit corresponding to ^ — > 00). The continuous 
and dashed lines display the corresponding numerical solution 
of PB theory. The dotted line is a guide for the eye indicating 
a slope Q = 1/2. The inset shows the same data on a linear 
scale. 



uation of the cell model. It should be kept in mind how- 
ever that l|5.8|l does not allow to define a Manning ra- 
dius for < ^ < ^Manningj whcrc the condensation phe- 
nomenon is already present, at least from a mathematical 
point of view with a change in the short distance behav- 
ior of y{r). In addition the scaling oc (o/k)^/^ valid 
for large enough ^ [see (|5.7ll ] is definition dependent. An 
alternative to (|5.8I) could be to define a characteristic ra- 
dius of the condensate through q{r* ) — ^c- Making use of 
(|4.13l) then implies that for a 1:1 salt nr* cx (ko)" with 
a = (arctan7r)/7r ~ 0.402 instead of a = 1/2 in (|5.7|) . As 
(|4.13(l . this is limited to high enough ^. For ^ = 1, the re- 
sults of section HV CI allow for an analytical computation 
of r* and yield a = {2/tt) arctan(7r/2) ~ 0.639. 

Our expressions also allow to discuss several quanti- 
ties that directly follow from the electric potential, such 
as the Bjerrum radius tb considered in Ref. [3^ and de- 
fined as the locus of an inflection point in the integrated 
counterion density when plotted as a function of radial 
distance. For 1:1 and 2:1 salts, this definition implies 
^(^b) = 1/2 while in the 1:2 case, we have g(r_B) = 1/4. 
Making use of Eqs. H5.5|l provide a transcendental equa- 
tion from which tb follows. 



B. Effective charges 

Of particular interest to describe interactions at large 
distances (typically r > as will be discussed in sec- 
tion lV C|l is the effective charge defined from the far field 
asymptotics (|1.3f) . This quantity is given by Eqs. H2.8|) 
where A follows from the expressions given in sections 
IIII Bl and unci The results pertaining to the limit 5, ^ 
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have been given in section HVAI but it is also possible to 
derive closed form relations in the saturation limit ^ oo 



'ssat 




10"^ 10"° 10"'' 10"* 10"' 10"' 10"' 



Kfl 

FIG. 8: Effective linear charge density at saturation, ^sat 
as a function of the radius 5 of the rod for a symmetric 1:1 
electrolyte. The full line is the analytical expression 15.9al l. 
the circles have been obtained by solving numerically Poisson- 
Boltzmann equation. The inset shows the same data on a 
linear scale. 

Figures |H1 and |51 show that the analytical expressions 
(|5.9|) are in good agreement with the numerical results for 
Ka — d < 0.1. For d > 0.1, deviations become apparent 
(see the inset of Fig. ^ . We will propose in section IVII 
an alternative approach that covers the whole range of 
Ka. 

We have also checked that the effective charges are cor- 
rectly described by the analytical predictions of section 
mil not only in the saturation regime but for arbitrary 
values of the charge ^ provided d < 0.1. 



C. Electric potential 

So far, we focused on several scalar quantities char- 
acterizing the electric potential, and thus the ionic dis- 
tribution. It is also instructive to compare the analyti- 
cal potentials to their "exact" PB counterparts. When 
^ > £,c, the relevant short scale expressions are those of 
Eqs. H3.7|l . We take the 2:1 situation as an illustrative 



where fi is given by H4.14|l . We therefore have (denoting 
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FIG. 9: Same as the main graph of Fig. |Hlfor 1:2 and 2:1 
electrolytes. The analytical expressions shown by the lines 
are those of Eqs. J5.9bt and Il5.9ct . The limiting values for 
K,a — > 0^ are shown by the dashed lines. 



example. As expected from the results of section IV Bl 
Figures [TUI and ITTI show that the far field behavior is well 
captured by H1.3|l . and it also appears that such an ex- 
pression may be used in practice for Kr > 1. On the 
other hand, expansion (|3.7c|) is very accurate at small 
distances, and may be extended up to nr of order 1. For 
the low values of kq such as that of Fig. 1101 expansion 
H^.Tcll produces the correct potential with an exceptional 
accuracy for several orders of magnitude in r. For higher 
values of ko, the short distance expansion (|X7c|) is still 
useful and reliable up to nr of order unity (see Fig. Ill|l . 
To illustrate the improvement over previous expansions, 
we come back to the 1:1 case in Fig. II 21 and compare our 
formula with the classic one of Ramanathan , recalled 
in H5.1|l and re-derived in the present work. Once the po- 
tential is known, one can compute the integrated charge 
q{r) = —{l/2)rdy/dr, displayed in Fig. ^1 The coin- 
cidence between the q versus logr inflection point and 
the Manning radius may be observed in the inset. The 
integrated charge -shown by the dotted line- following 
from H5.1|l cannot reproduce such a feature, and we also 
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FIG. 10: Electrostatic potential in a 2:1 electrolyte as a func- 
tion of radial distance for a rod with a high bare charge 
(5 ~ 11) and Kfl = 10~^. The circles represent the numerical 
solution of PB theory, the dotted line is for the short dis- 
tance formula 13.7cl l and the dashed line is for the far field 
expression 11.31 . In the inset, the same results are shown on 
a Unear-log scale. 




FIG. 11: Same as Fig. [TOlwith 
linear scale. 
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FIG. 12: Comparison of the numerical solution of PB theory 
with expansion Il3.7al shown by the dotted line and with the 
formula of Ramanathan [see 115.111 and dashed line] for nia — 
10"^, ^^25 and a 1:1 sah. 




FIG. 13: Integrated charge as a function of rescaled distance 
for the same parameters as in Fig. 1121 (symbols and lines 
have the same meaning here) . The inset shows the same data 
in the vicinity of the predicted Manning radius, indicated by 
the arrow. The prediction follows from l|5.6a^ or equivalently 
15.7all since we consider here a highly charged rod. 



observe that the most significant decay of q(r) from its 
initial value as requested by Gauss theorem, i.e. 25 in 
Figure [T5)) to unity occurs on a much shorter scale than 
the Manning radius. 



VI. EFFECTIVE CHARGE AT SATURATION 
AND THE LARGEST EIGENVALUE OF THE 
OPERATORS Ka 

The results of the previous sections are broadly speak- 
ing restricted to the regime ko < 1. For the specific 
problem of the effective charge in the 1:1 situation, it is 
however possible to obtain expressions that hold for all 
salinities. 



The effective charge at saturation, for any arbitrary 
value of a, has an interesting relation with the largest 
eigenvalue of the operators defined in Eqs. (|2.3|l and (|2.6|l . 
To see this, consider first the 1:1 electrolyte, and the 
solution H2.2|l to Poisson-Boltzmann equation. As the 
bare charge ^ increases the parameter A (related to the 
effective charge ^eff) increases. At saturation (^ +oo), 
A = Asat is such that the electric field obtained from the 
solution (|2.2|l diverges at f = a. Clearly this happens if 
det(l — XKa) = 0, thus if A is the inverse of an eigenvalue 
of ii'aE3- Now, since when ^ = 0, A = and it increases 
as ^ increases, it appears that at saturation A = Asat is 
equal to the inverse of the largest eigenvalue of Ka ■ 

The same analysis applies to the 1:2 and 2:1 elec- 
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trolytes. For a 1:2 electrolyte, for positive saturation 
^ ^ +00, Agat is the inverse of the largest eigenvalue of 
, whereas for a 2:1 electrolyte, when ^ +00, Asat 
is the inverse of the largest eigenvalue of —K^\ 

We have not been able to find explicitly the eigenvalues 
of the operators Ka, K^'^ and k''^\ for any arbitrary 
value of a. However one can use approximate methods 
to find estimates for the largest eigenvalue. 

For the 1:1 electrolyte, it is shown in the 
appendix of [s^ . that for any function G 
L^(0, 00, e^"^"+" ^/^d-u), the largest eigenvalue Ag^J: of 
Ka satisfies 



< A: 



(6.1) 



where (• , •) is the scalar product of 
L2(0,cx),e-°(''+""')/2^^)^ Using the test function 
(j){u) = 1/ \/u, we obtain 



A« 



< 



7rr(o,2a) 



(6.2) 



with Kq the modified Bessel function and r(0, z) = 

+00 



*/tdt. Using Eq. H2.8a(l . this finally gives an up- 



per bound for the effective charge at saturation, for a 1:1 
electrolyte, for any arbitrary value of a, < "p, 
with 



sat, up 



2aKi{d)Ko{d) 
7rr(0,2a) 



(6.3) 



As a ^ 0, 2/7r, thus has the same limit as the 

exact £_lg. This is expected since in the appendix [s^ . 
this upper bound was used to prove that the supremum 
of the largest eigenvalue of Ka as a ^ is tt. 
Interestingly, when a 3> 1 we have 



sat, up 



25 



(6.4) 



which is the same asymptotic behavior of the true effec- 
tive saturated charge obtained in [s^ by a different 
approach than the present one. It is the effective charge 
of an infinite plane (2Ka) plus a correction (3/2) that is 
obtained using a small curvature expansion. 

Figure shows a comparison between "p from 
Eq. 1)6. 3|l and the effective charge at saturation ob- 
tained numerically. Surprisingly, it turns out S,^^*"' "p is 
not only an upper bound for but also a very good 
estimate for for any value of a. If a > 10^^ it is in 
very good agreement with the numerical data. However, 
for very small radius a <C 1, the estimate (|5.9a|l of sec- 
tion is better than the ansatz (|6.3|) (see the inset of 
Fig. 1141 and compare it to Fig. |SJ). 



Vn. CONCLUSION 

The mathematical results derived in the framework of 
Painleve/Toda type equations provides much insight into 
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FIG. 14: The upper bound 16.311 used as an ansatz for the 
effective charge at saturation as a function of the radius a of 
the rod for a symmetric 1:1 electrolyte. The full line is the an- 
alytical expression 16.31 , while the circles have been obtained 
by numerical resolution of Poisson-Boltzmann equation. 



the screening behavior of microions in the vicinity of a 
charged rod-like macroion. Although the mapping be- 
tween the Poisson-Boltzmann equation for a certain class 
of electrolytes (1:1, 1:2 and 2:1) onto a Painleve/Toda 
equation is in itself not new, previous analysis were con- 
cerned with the no salt limit where the ratio of macroion 
radius a to Debye length k^^ vanishes, and the practical 
consequences and implications of a finite salt concentra- 
tion had not been drawned. We have shown here that 
systematic logarithmic, and thus strong, corrections arise 
for various quantities of interest. All finite salt results 
reported here are new, and significantly improve previ- 
ously available expressions. In addition, the 2:1 situation 
worked out here had not been considered before, even in 
the limit of vanishing salt content. 

The experimental relevance of our findings has been 
sketched in the introduction. A rough criterion for 
the validity of our mean-field approach is that F = 
z^l"^ ^ ^Ib I i^'na) < 2 where z denotes the valency of 
counter-ions. For 1:1 and 2:1 electrolytes, the criterion 
is fulfilled even by the most highly charged polymers, 
such as double-stranded DNA for which a ~ 1 nm and 
^ ~ 4, at least in water where £b — 0.7 nm. In the 1:2 
case, we obtain for dsDNA F ~ 2, for which it is diffi- 
cult to weight the importance of microionic correlation 
against the mean-field effects reported here. In any case, 
our results apply to single stranded DNA and polymers 
of lower line charge such as hyaluronan |39l |. An upper 
limit should also be put to salt content {k£b < 1), but 
such a requirement is irrelevant for our purposes since 
only the low salt regime has been discussed. 

Section IIIII contains the core of our analysis. It ap- 
peared there that the short scale behavior of the electric 
potential depends on the position of the reduced bare 
charge ^ with respect to a threshold ^c- One may con- 
sider that this threshold is a remnant of its vanishing 
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salt counterpart, which signals the onset of counter- ion 
condensation. It should be kept in mind however that 
counter-ion condensation is not an all or nothing process, 
but a gradually built phenomenon. 

The charge is strongly salt dependent, and coincides 
with Manning critical value in the limit — > 0. On the 
other hand, the large distance behavior is always of the 
same functional form [see H1.3|l ] and is governed by the 
effective charge ^eff- This latter quantity depends on a 
parameter A that also plays a crucial role in the descrip- 
tion of the short distance behavior. We have proposed a 
simple approximation to compute explicitly A for 1:1, 2:1 
and 1:2 electrolytes, from which all other quantities fol- 
low. The analytical predictions have been tested against 
the numerical solution of Poisson-Boltzmann theory, and 
shown to be remarkably accurate [of particular practical 
interest for highly charged polyelectrolytes are Eas. H3.7|) 
and l|1.3(l with n given by l|4.13|l and where A is given as a 
function of ^ in section Fill Bj . In particular, simple ana- 
lytical results have been derived for the Manning radius, 
that is often used to quantify the lateral extension of the 
condensate that may form around a cylindrical polyion. 
A few other measures of the condensate thickness have 
been commented upon. 

Our approach is free of the matching procedures 
|40l I41L or ad-hoc though educated assumptions 
underlying previous work, and therefore provides expan- 
sions with controlled error for the electric potential. The 
analysis of sections IIIII to Ivl requires that a < k~^. For 
larger salinities, contributions that have been neglected 
here become relevant. However, the spectral analysis of 
section fVll provides for all values of Ka an upper bound 
for the far field signature of the polyion, from which an 
excellent approximation of the effective charge at satura- 
tion may be obtained. 

The analytical results obtained emphasize the essential 
difference between the critical Manning charge for con- 
densation and effective charge of the macroion, the latter 
being defined from the large scale electric behavior. A 
consequence is that within a simplified two state model 
where the population of microion is divided into a con- 
densed region and a free population, free ions cannot be 
treated within a linearized theory, at variance with com- 
mon belief and practice. Consistency with exact results 
requires that non linear effects are still at work in the 
free region and significantly decrease the effective charge 
compared to the critical one. In the 1:2 case at ultra low 
salt, for all bare charges ^ larger than ^c, the ratio of 
Coff to (with equal to ^Manning bere since na 0) 
is equal to a/S/tt ~ 0.55. This ratio is closer to unity 
for 2:1 electrolytes [3-\/3/(27r) ~ 0.83], and intermediate 
for a 1:1 salt {2/tt ~ 0.64). We therefore emphasize that 
the distinction between bound and free ions is not only 
arbitrary but also potentially misleading. 
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APPENDIX A: PAINLEVE CLASSIFICATION, A 
BRIEF REMINDER 

Polynomial non-linear differential equations of the 
form 

A(.,,)g-^B(.,,)| + C(.,,)(g)%i.(.,,) = 0, 

(Al) 

where the functions A,B,C,D are polynomial in y and 
analytic in x, have been classified with respect to the 
character of the singular points of the solutions. Of spe- 
cial interest are the equations for which branch points 
and essential singularities do not depend on initial con- 
ditions (hence the only movable singularities are poles). 
Fifty canonical types of equations with the above prop- 
erty have been uncovered, most of which (44) are in- 
tegrable in terms of elementary functions. Solving 
the remaining 6 types requires the introduction of new 
(Painleve) transcendental functions. The third mem- 
ber (Painleve III) of this family of 6 corresponds to the 
generic form 

(Py f dy\^ dy , 4 / a 

xy—r = a; — - y— + ax + by + cy +dxy . (A2) 
dx^ \dx J dx 

APPENDIX B: SHORT DISTANCE BEHAVIOR 

FOR A > Ac 

For a highly charged rod, when ^ > ^c, we have 
5Re(A) — ^Manning- Then, the last "higher" order term in 
Eqs. (|3.5|l become of the same order as the 0(1) terms. 
Thus the small-f asymptotics of y{f) are different, when 
^ > t han the ones given by Eqs. H3.5|l . 

In |22| the asymptotics for A > 1 in the 1:1 case were 
studied and in (3| the ones for the 1:2 case for A > Ac 
were obtained. We will not reproduce those calculations 
here, but to illustrate the method from [2^ W^. we will 
compute the small-r asymptotics for the 2:1 case, for A > 
Ac , which has not been previously considered. 

The asymptotic form (|3.5c|) is actually valid even if A 
is complex, provided that it satisfies (|3.14|) . To study 
the asymptotics when A > Ac we can consider that A 
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approaches the cut [Ac, +00) from below, for instance. 
Then we can rewrite (j3.13b|l as 



with 



In 



— cosh 



A = I in 

2 ^ 



3A 1 
2X'r~ 2 



3A 

2a: 



-1 f 

2Xr 



(Bl) 



(B2a) 



(B2b) 



Replacing into (|3.5c|) and keeping only the first two dom- 
inant terms (which arc of the same order) gives 



-V21/2 



with 



/ r 



V6\/3 



r(i-»f)r(i-zM) 
r(i + if)r(i + i/i) 



1/2 



(B3) 



(B4) 



If one chooses A to approach the cut [Ac, +00) from above 
then /i is changed into and the final result is un- 
changed. 

Finally, the potential 2/21 is given by 



,-2,21/2 _ 



2f 
3/iV6 



— ;r In — 1= 
2 6%/3 



where 



= <{ in[r(i - ^)r(i - z^)] 



(B5) 



(B6) 



For the sake of completeness, we reproduce here small- 
r asymptotics for the 1:1 and the 1:2 cases which were 
computed in Refs. (23 and respectively. For the 

1:1 electrolyte, when A > Ac'^ '^'', the electric potential is 
given by 



sin 

4// 



-2^1n^ -|-2*(i^i)(^) 



(B7a) 



with 



and 



At = - cosh"'(7rA) > (1:1) (B7b) 

TT 



*(i^i)(Ai) =5m[lnr(l + vi)] 



(B7c) 



For the 1:2 electrolyte, when A > Ac ' , the asymptotics 
are Hi 



g-yi2(f) ^ ^sin 



-3/iln^-2vI/(i^2)(^) 



now with 



1 _i /I A , 



and 

=Qto<| In 



r ( ^-^ 1 r(l - z^i) 



+0(f^) (1 
(B8) 

(B9) 



(BIO) 



For practical purposes, since /i is at most of order 
1/1 lna|, we can expand the functions 



^,[--^ + \n2) + o{^x') 



(Bll) 
(B12) 
(B13) 



The asymptotics presented in Eqs. I|3.7(l in the main text 
use this approximation for ^{^). 

In principle this means that our expressions in the main 
text for jjL are accurate only to order 0(1/1 In a|^). How- 
ever, an explicit computation using further terms in the 
Taylor expansion of \E'(/x) shows that our results, in the 
form presented in the main text, are actually accurate 
up to order l/|lna|^ for ^ and order l/|lna|'* for the 
effective charges. 

To illustrate this, consider, for example, the value 
of /Zsat for the 1:1 case. Using a Taylor expansion of 
\I''^'^'(At) up to order gives 



(1:1) 
Msat 



21n§ 



7 
In I 



7' I'-^j^^Hi) 



(1^4) 



0((lna)- 



(B14) 



where '4'{x) = d\nT{x)/dx is the digamma function and ijj'^'^^x) its second derivative. Replacing this expression into 
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Eqs. (|3.12() and H2.8|l . yields the effective charge at saturation 



7^ 4(ln|)^ 2(ln|)^ 192 (in f)"* 



(1-1)^ 



o((ina) 



(B15) 



A direct comparison of Eas. (|B15|) and (|B14|) with H4.14a|) presented in the main text are indeed accurate up to 
and H5.9a|l shows that the (more compact) expressions order l/|lnap for ^ and order l/|lna|'* for ^sat- 



A. Khokhlov in "Soft and Fragile Matter" , edited by M.E. 
Gates and M.R. Evans, Scottish Graduate Textbook se- 
ries, loP (2000). 

G.S. Manning, J. Ghem. Phys. 51, 924 (1969). 

F. Oosawa, Polyelectrolytes, Dekker, New York (1971). 
T. Alfi-ey Jr, P.W. Berg and H. Morawetz, J Polym. Sci. 
7, 543 (1951); R.M. Fuoss, A. Katchalsky and S.F. Lifson, 
P. Natl. Acad. Sci. USA 37, 579 (1951). 

M. Fixman, J. Ghem. Phys. 70, 4995 (1979). 

M. Gueron and G. Weisbuch, Biopolymers 19, 353 

(1980). 

G. V. Ramanathan, J. Ghem. Phys. 78, 3223 (1983). 

M. Le Bret and B.H. Zimm, Biopolymers 23, 287 (1984). 
J.S. McGaskill and E.D. Fackerell, J. Ghem. Soc, Fara- 
day Trans. 2 84, 161 (1988). 

A. L. Kholodenko and A.L. Beyerlein, Phys. Rev. Lett. 
74, 4679 (1995). 

G. A. Tracy and H. Widom, Physica A 244, 402 (1997). 
Y. Levin, Physica A 257, 408 (1998). 

M. Deserno, G. Holm and S. May, Macromolecules 33, 
199 (2000). 

P.L. Hansen, R. Podgornik and V.A. Parsegian, Phys. 
Rev. E 64, 021907 (2001). 

B. O'Shaughnessy and Q. Yang, Phys. Rev. Lett. 94, 
048302 (2005). 

L. Belloni, J. Phys.: Condens. Matter 12, R549 (2000). 
J. -P. Hansen and H. Lowen, Annu. Rev. Phys. Ghem. 51, 
209 (2000). 

A. Y. Grosberg, T.T. Nguyen and B.I. Shklovskii, Rev. 
Mod. Phys. 74, 329 (2002). 

Y. Levin, Rep. Prog. Phys. 65, 1577 (2002). 

R.R. Netz, Eur. Phys. J. E 5, 557 (2001); A. Naji and 

R.R. Netz, Eur. Phys. J. E 13, 43 (2003). 

Y. Levin, E. Trizac, and L. Bocquet, J. Phys.: Condens. 

Matter 15, S3523 (2003). 

B. M. McGoy, G.A. Tracy and T.T. Wu, J. Math. Phys. 
18, 1058 (1977). 

H. Widom, Gommun. Math. Phys. 184, 653 (1997). 
G.A. Tracy and H. Widom, Gommun. Math. Phys. 190, 
697 (1998). 

E. Trizac and G. Tellez, Phys. Rev. Lett 96, 038302 
(2006). 

see e.g. H.T. Davis, Introduction to nonlinear differential 

and integral equations, Dover 1962. 

G.J. Benham, J. Ghem. Phys. 79, 1969 (1983). 

L. Samaj and I. Travenec, J. Stat. Phys. 101, 713 (2000). 



[29] F. Gornu and B. Jancovici, J. Stat. Phys. 49,33 (1987). 

[30] T. Odijk, Ghem. Phys. Letters 100, 145 (1983). 

[31] E. Trizac, L. Bocquet, M. Aubouy and H.H. von 

Griinberg, Langmuir 19, 4027 (2003). 
[32] G. Tellez and E. Trizac, Phys. Rev. E 70, 011404 (2004). 
[33] G. Tellez and E. Trizac, Phys. Rev. E 68, 061401 (2003). 
[34] R.R. Netz and J.-F. Joanny, Macromolecules 31, 5123 

(1998). 

[35] G. Weisbuch and M. Gueron, J. Phys. Ghem. 85. 517 
(1981). 

[36] H. Qian and J.A. Schellman, J. Phys. Ghem. B 104, 
11528 (2000). 

[37] G.A. Tracy, Gomm. Math. Phys. 142, 297 (1991). 

[38] M. Aubouy, E. Trizac and L. Bocquet, J. Phys. A: Math. 

Gen. 36, 5835 (2003). 
[39] E. Buhler and F. Boue, Macromolecules 37, 1600 (2004). 
[40] A.D. MacGillivray and J.J. Winkleman Jr, J. Ghem. 

Phys. 45, 2184 (1965). 
[41] J.R. Philip and R.A. Wooding, J. Ghem. Phys. 52, 953 

(1970). 

[42] H. van Keulen and J. A.M. Smit, J. GoUoid Interface Sci- 
ence 170, 134 (1995). 

[43] R. Tuinier, J. GoU. Interface Sci. 258, 45 (2003). 

[44] Note that the 1 : 2 case with ^ > 0, with colons of valency 
z\ = 4-1 and counterions of valency 22 = —2, is identical 
to the case z\ = +2, Z2 = —1 with ^ < 0. Equation II. let 
and the boundary condition 11.21 1 is of course formally 
equivalent to Eq. ll.lbl after the substitutions y — > —y 
and ^ 

[45] In this respect, there is a further approximation in 113.1711 
compared to 13.1611 since it relies on 13.7cL which is not 
exact. A similar remark holds for 13.191 or 13.2111 . It is 
however a simple task to get rid of this extra degree of 
approximation using the expressions given in Appendix 
ini In any case, the resulting relation between ^ and fj, 
is only asymptotically exact in the limit a —f 0, and re- 
sorting to the full expressions of Appendix^brings very 
little improvement in terms of numerical accuracy. 

[46] We also note that references |9, 10] suffer from several 
typographic errors that make detailed comparison quite 
difficult. 

[47] This analysis is for (positive) saturation when ^ > 0. In 
the case ^ < 0, the saturation is reached when det(l -|- 
XKa) = 0. 



